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AEROSOL  SAMPLE  ACQUISITION 

FOR  CHEMICAL  AND  BIOLOGICAL  AGENT  DETECTION 

i.  Introduction 

The  objective  of  this  project  is  to  investigate  means  for 
sampling  aerosol  particles  from  the  ambient  atmosphere  for  chemical 
and  biological  detection.  Sample  acquisition  for  such  purpose  require 
innovative  approaches  involving  sample  extraction  from  the  ambient 
atmosphere,  means  for  concentrating  the  aerosol  particles  from  the 
extracted  sample  into  a  smaller  airstream,  collecting  the  concentrated 
particles  onto  a  suitable  substrate  and  then  volatilizing  and/or 
pyrolyzing  the  collected  sample  for  subsequent  physical  or  chemical 
analysis.  While  considerable  knowledge  is  available  on  sample 
extraction  from  the  ambient  atmosphere,  less  is  known  about  sample 
concentration  and  volatilization  and  pyrolyzing  methods.  The  major 
focus  of  this  project  is  to  perform  theoretical  and  experimental  study 
of  sample  concentration  and  volatilization/pyrolysis  methods  with 
the  goal  of  developing  sufficient  understanding  of  the  processes  so 
that  complete  sample  concentration  and  volatilization/pyrolysis 
systems  can  be  developed  to  suit  different  needs. 

The  report  is  divided  into  several  sections  with  sections  2  -  4 
devoted  to  the  theoretical  and  experimental  study  of  virtual 
impactor  concentrators.  Following  a  review  of  the  theoretical 
methods  used  for  virtual  impactor  studies  in  Section  2,  a  design 
analysis  of  virtual  impactor  concentrator  is  made  in  Section  3.  A 
specific  virtual  impactor  concentrator  with  a  design  flow  rate  of  22 
liters  per  minute  has  been  designed,  constructed  and  tested  and  this 
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is  described  in  detail  in  Section  4.  Results  of  this  study  provides 
information  for  virtual  impactor  concentrator  design  of  nearly  any 
flowrate  and  particle  cut  size  in  the  future.  In  Section  5,  the 
problems  related  to  particle  collection  for  volatilization  and  pyrolysis 
are  studied  and  analyzed.  Several  experimental  volatilizer/pyrolyzer 
concepts  were  analysis  and  compared.  Finally,  an  experimental 
volatilizer/  pyvolyzer  was  constructed  and  tested.  On  the  basis  of  the 
study,  a  new  volatilizer/pyrolyzer  was  then  proposed  for  possible 
future  development  for  various  applications  as  described  in  Section 
6. 

2.  Theoretical.  fia&kgiQund.,  and  Analysis 

One  of  the  most  powerful  techniques  to  determine  the  details 
of  fluid  flow  Fields  and  particle  trajectories  in  an  aerosol  system  is  by 
use  of  numerical  methods  to  solve  the  equations  governing  the  flow 
field  and  the  particle  motion  in  that  flow  field.  The  governing 
equations  for  the  flow  field  are  the  Navier-Stokes  equations,  and 
those  for  the  particle  trajectories,  the  equations  of  motion  of  the 
particles  derived  from  the  Newton's  law  relating  the  forces  acting  on 
the  particle  and  the  particle  inertia. 

The  first  part  of  this  technique  (flow  field  analysis)  has  been 
used  by  many  investigators  to  determining  the  details  of  the  flow 
around  objects  or  through  confined  spaces.  However,  we  have  used 
both  parts  of  the  technique  extensively  to  determine  the  particle 
trajectories  as  well  as  the  flow  field,  and  in  some  cases  have  used  the 
technique  to  determine  the  particle  concentration  profile  in  an 
aerosol  system.  The  technique  has  been  used  in  detailed  studies  of 
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particle  motion  in  inertial  impactors,  filters,  elutriators,  clean  rooms, 
glove  boxes,  at  inlets  as  well  as  in  virtual  impactors.  In  many  of 
these  studies  the  theoretical  results  have  been  compared  with 
experimental  results  and  agreement  has  been  found  to  be  excellent  if 
the  boundary  conditions  used  in  the  actual  system  closely 
approximates  those  used  in  the  theoretical  analysis. 

2.1  The  virtual ,  impactor 

In  a  conventional  impactor  a  jet  of  particle  laden  air  is 
accelerated  through  a  nczzle  and  directed  at  a  flat  plate.  Particles 
larger  than  a  certain  critical  size,  the  so-called  cut  size  of  the 
impactor,  would  impact  upon  the  plate  while  smaller  particles  would 
follow  the  airstream  and  escape  impaction.  In  a  "virtual  impactor"  an 
accelerating  nozzle  is  used  as  in  the  conventional  real  impactor  to 
accelerate  the  airstream  to  a  suitably  high  velocity,  but  the  impaction 
plate  is  replaced  by  a  void  space  in  the  form  of  a  receiving  tube 
through  which  a  small  fraction  of  the  flow  is  allowed  to  pass.  The 
main  flow  is  deflected  by  the  receiving  tube,  but  the  large  particles 
are  able  to  cross  the  streamlines  and  enter  the  receiving  tube  on 
account  of  their  inertia,  while  the  small  particles  are  carried  by  the 
deflected  airstream  around  the  impactor.  The  output  of  the  virtual 
impactor  are  in  the  form  of  two  airstreams.  The  small  stream  (the 
minor  flow)  contains  the  large  particles  inertially  impacted  into  the 
void  space  of  the  receiving  tube,  and  the  main  airstream  (the  major 
flow),  the  small  particles  contained  in  the  original  flow.  In  contrast 
to  the  conventional  real  impactor  where  the  large  particles  are 
impacted  on  a  solid  surface,  the  virtual  impactor  has  the  large 
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particles  impacted  into  a  void  space  and  these  particles  would  thus 
remain  suspended  or  airborne. 

The  virtual  impactor  has  two  features  which  are  desirable  for 
the  continuous  sampling  and  analysis  of  airborne  particles.  First, 
both  the  large  and  small  particles  remain  airborne  and  can  be  easily 
transported  to  another  instrument  for  collection  (such  as  a  real 
impactor),  to  an  instrument  for  analyses  (  such  as  an  optical  counter 
for  continuous  size  distribiuion  monitoring)  or  be  removed  from  the 
airstream  by  filtration.  Second,  the  particles  larger  than  the  cut  size 
of  the  impactor  are  concentrated  into  the  small  fraction  of  the  flow 
passing  through  the  receiving  tube.  This  makes  analysis  easier  and 
more  sensitive  than  if  the  particles  were  in  the  total  flow  entering 
the  impactor. 

The  major  features  of  the  flow  in  a  virtual  impactor  is  rather 
simple.  The  air  exits  the  nozzle  and  penetrates  into  the  receiving 
tube  that  is  about  30%  larger  in  diameter  than  the  nozzle.  The  depth 
(measured  in  number  of  receiving  tube  diameters)  to  which  the  air 
penetrates  into  the  receiving  tube  depends  upon  the  Reynolds 
number  of  the  air  passing  through  the  nozzle.  The  Reynolds  number 
is  based  on  the  nozzle  diameter  and  the  average  velocity  of  the  air  in 
the  nozzle.  Most  of  the  air  reverses  direction  inside  the  receiving 
tube  and  exits  the  receiving  tube  at  its  outer  diameter.  Particles 
with  large  enough  inertia  cannot  reverse  direction  and  penetrate 
farther  into  the  receiving  tube.  Since  a  small  fraction  of  the  flow 
entering  the  receiving  tube  is  drawn  through  the  tube,  the  large 
particles  are  concentrated  into  this  small  flow. 
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Although  the  flow  in  the  virtual  impactor  is  rather  simple  in  its 
gross  features,  there  is  no  analytical  solution  for  the  basic  flow  field 
and  particle  trajectory  equations.  Therefore,  the  problem  must  be 
solved  numerically  to  determine  the  details  of  the  flow  field  and 
particle  trajectories. 

2.2  Numerical  Techniques  for  Flow  Field  Analysis 

There  are  several  numerical  routines  for  solving  the  Navier- 
Stokes  equations.  One  routine  expresses  the  Navier-Stokes  equations 
and  the  continuity  equation  in  terms  of  the  the  stream  function  and 
the  vorticity.  These  two  quantities  are  then  solved  numerically  at  the 
node  points  covering  the  computation  domain.  The  velocity 
components  at  the  node  points  are  then  detrrmined  from  knowledge 
of  the  value  of  the  stream  function  at  these  points.  Other  routines, 
such  as  SIMPLE  and  SIMPLER,  solves  directly  for  the  velocity 
components  and  pressure  at  the  node  points.  Both  types  of  routines 
should  give  the  same  answer  since  the  same  basic  sets  of  equations 
derived  from  the  same  basic  Navier-Stokes  equations  are  solved.  In 
either  case,  the  values  of  the  velocity  components  at  the  node  points 
is  the  final  result.  By  using  these  velocity  components,  the  particle 
trajectories  are  then  calculated  in  the  same  manner  for  both  cases. 

We  have  used  both  types  of  routines.  However,  because  the 
boundary  conditions  are  slightly  easier  to  define  for  the  stream 
function  -  vorticity  routine,  we  have  chosen  it  to  be  the  routine  for 
analysis  of  the  virtual  impactor. 

The  general  method  of  solution  is  to  first  express  the  Navier- 
Stokes  equations  in  terms  of  the  vorticity  and  the  stream  function. 
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The  resulting  differential  equations  are  then  changed  to  a  finite 
difference  form  and  solved  by  the  method  of  relaxation  over  a  grid 
of  node  points  covering  the  field  of  interest.  This  method  of  analysis 
follows  closely  the  general  solution  procedure  described  by  GOSMAN, 
et  al.  (1969). 

In  the  method  devised  by  GOSMAN,  et  al.  (1969)  it  is  not 
necessary  to  have  equally  spaced  grid  lines.  Therefore,  more  node 
points  can  be  placed  in  areas  where  the  velocity  gradients  are  the 
largest.  It  is  also  assumed  that  the  density  and  viscosity  are  constant 
and  that  there  *  o  swirl  so  that  there  is  no  0  component  of  velocity. 
Considering  these  assumptions  the  Navier-Stokes  equations  are: 

r  direction: 


,  avr'  av;  3p  32Vr  !3V'r  V'r  a2yr 


z  direction: 


z  av*  av* 


§2.  /d2v*z  1 5V'Z  52v* 

=  ‘az’  +  ^ar'2  +r*  3r*  + 


(2) 


and  the  continuity  equation  is: 


av'r  V*r  av*z  _ 
ar*  +  r*  +  9z*  “° 


(3) 


m causes *sbs j.ts 
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The  continuity  equation  is  satisfied  by  the  stream  function,  y’, 
defined  as  follows: 


_  L*jL  v.  _!^ 

“V  dr”  r  ~r’  dz' 


(4) 


The  vorticity,  co\  is  defined  as: 


O) 


avz  avy 

dr'  ‘  dz' 


(5) 


The  above  five  equations  are  combined  into  two  differential 
equations  which  fully  describe  the  flow  field  in  terms  of  to’/r’  and  the 
stream  function.  These  equations  are  made  dimensionless,  as 
denoted  by  dropping  the  primes,  by  means  of  the  jet  diameter,  W, 
and  the  average  air  velocity  at  this  diameter,  V0.  The  dimensionless 

terms  are: 


•1 

II 

3h 

N 

II 

3  b 

(6) 

¥ 

V  =  V0W2 

(7) 

co  ta^rl  - 

(8) 

r  =v0/w2 

and 
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Re  = 


pv0w 

4 


(9) 


The  Navier-Stokes  equations  (equations  1  and  2),  continuity 
equation  (equation  3)  and  the  definition  of  the  stream  function 
(equation  4)  and  the  vorticity  (equation  5)  have  been  combined  into 
two  differential  equations  which  fully  describe  the  flow  field  in 
terms  of  w/r  and  the  stream  function.  These  two  equations  are 
made  dimensionless  by  means  of  equations  (6)  to  (9)  and  the  results 
are: 


a( 

r2- 


£0  d_Vj/x 

r  dz' 
dr 


*7 


“IS)  d(r3 


dr' 


-  r“ 


dz 


-hi- 


dto/ts 

dr  ' 


*3ir) 


dr 


dz 


■} 


(10) 


and 


a  _  9/1  9\|/x  d_  ,  1  dv/s 
r  r  “  8r  ^r  9r^  dz  ^r  dz ' 


(ID 


Equations  (10)  and  (11)  are  integrated  over  a  finite  area 
around  a  typical  node  point  and  two  finite  difference  equations  are 
obtained.  The  typical  node  point  is  denoted  by  the  subscript  P  and 
the  surrounding  node  points  by  E,  W,  N,  S,  SE,  NE,  SW  and  NW. 
Figure  1  shows  the  orientation  of  these  node  points  with  respect  to 
the  r  and  z  directions. 
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The  convective  terms  (the  2  terms  on  the  left-hand  side  of 
equation  (10)  are  integrated  using  the  upwind  difference  method. 

The  result  of  integrating  equation  (10)  and  solving  for  ca/r  at 
point  P  is: 


(co/r)E(AE  +  RBE)  +  (to/r)w(AW  +  RBW)  +  (co/r)N(AN  + 
(“)p  =  fp(RBN)  +  (o)/r)s  AS  +  rp  (RBS) 


AE  +  AW  +  AN  +  AS  +  RBE  +  RBW  +  r2p  (RBN  +  RBS) 

(12) 


where: 

AE  =  {(vS  +  vSE  -  vN  -  vNE)  +  lyS  +  ySE  -  yN  -  yNE !}  DV 
AW=  {(yN  +  yNW  -  yS  -  ySW)  +  lyN  +  yNW  -  yS  -  ySW  1}  DV 
AN  =  {(yE  +  yNE  -  yW  -  yNW)  +  lyN  +  yNW  -  yS  -  ySW  1}  DV 
AS  =  {(yW  +  ySW  -  yE  *  ySE)  +  lyW  +  ySW  *  yE  -  ySE !}  DV 


(13) 


DV  = 


_ Is _ 

2(rE  *  rw)(2N*zS) 
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RBE 


.ill  _ 1 _ 1  3 

4  Re  rp  (rg  -  rp  )  (rg  -  ry/  )  ^  +  V 


RBW=i  ib  i  (TF^  (rw  +  rp)3 


RBN  =  J _ l _ I _ 

^Re  (z^  -  zp)  (zjj  •  z§) 


(14) 


BM  „  1  1 _ 1  . 

RBS  —  2p  /  _  \  _  .  ) 

Ke  (Zp  -  ZS)  zjsj  -  Zg 


The  result  of  integrating  equation  (11)  and  solving  for  the 
stream  function  at  point  P  is: 


\|/E  BE  +  \|/\y  BW  +  BN  BS  +  (to/r)p 

w  - - ^2 - XJ> - 

VP  “  __  1 


BE  +  BW  +  — 2  (BN  +  BS) 
rP 


(15) 


where: 


nr?  ,  1  1 
BE  =  4  ~ 


rp  (rE  *  rp)  (rE  •  rW  )  (rE  +  rp) 


BW  =  4  —  1 


rp  (rp  -  rW  )  (rE  *  rW  >  (rp  +  rW ) 


(16) 
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BN  =  2 ; - r  7 - r 

(ZN  *  zp>  (ZN  ‘  ZS> 

1 _ 1 

(zp  -  zs)  (zjsj  -  zg) 


Equations  (12)  and  (15)  are  the  equations  used  in  calculating 
co/r  and  y  at  each  node  point. 

Iteration  equations  are  also  required  for  grid  points  on  the 
boundary.  These  equations  for  the  boundary  node  points  are 
somewhat  different  than  for  the  interior  nodes.  For  each  of  these 
boundary  nodes  there  is  a  condition  for  the  stream  function  and  the 
vorticity. 

First,  the  stream  functions  along  the  surface  defining  the  nozzle 

is  held  constant  and  made  equal  to  zero.  Also,  along  the  centerline 
the  stream  function  is  held  constant  at  Vmax>  where  Vmax  =  1/8. 

Finally,  the  stream  function  at  the  receiving  tubes  wall  is  constant 
and  equal  to  the  percentage  of  Vmax  that  is  equal  to  the  fraction  of 

the  flow  in  the  major  flow. 

At  the  entrance  Vz  is  assumed  to  be  constant  so  that; 


.£  t  W 

V  ¥max  2  '  ^entrance 


(17) 


The  boundary  condition  at  all  exits  assume  the  streamlines  are 
straight  or: 


¥IN  ~  ¥INM 


(18) 
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Where  subscript  IN  denotes  boundary  nodes  and  INM  denotes  the 
adjacent  interior  node. 

The  vorticity  boundary  conditions  along  a  solid  surface  can  be 
deduced  from  the  fact  that  the  velocities  along  the  wall  and 
perpendicular  to  the  wall  are  set  equal  to  zero  and  the  gradients 
perpendicular  to  the  wall  are  assumed  to  be  much  larger  than  those 
along  the  wall.  Therefore,  for  surfaces  of  constant  z,  equation  (10) 
becomes: 


02(co/r)  _ 
dz 2 


(19) 


and  the  equation  (11)  becomes: 

o)  1  32v 
r  r2  9z2 


(20) 


leading  to  the  boundary  condition: 


3  ^NP  "  (®/r)jqp 

rp  (zp  -  ZNP)2  2 


(21) 


where  the  subscript  P  denotes  any  node  on  the  surface  and  the 
subscript  NP  denotes  the  adjacent  interior  node. 

For  surface  where  r=constant,  equation  (10)  becomes: 
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equation  (11)  becomes: 

co  _  .  1_  B_  ,1 
r  ‘  r  3r  '  r 


(22) 


(23) 


and  the  boundary  condition  becomes: 


3  4 

/OK  ,GK  rp  1  /rNP\  rNP 

(r)p  =  (VNP  -  Vp)  +  (r)Np  (rp  -  rNp)  1  6  [  4  rNP  ln  (  rp  ^  ‘  Ip2  + 

- ^ - L  r4  2  ln  . 1 +  r  21  -  -  (r  2  -  r  2)2 

P’  (rp  -  rNP)  1  6  1  rNP  ln  V  rp  /  2  +  rP  J  8  (rP  rNP  } 


(24) 


Where  the  subscripts  P  denote  grid  points  on  a  surface  and  NP 
denotes  the  adjacent  interior  grid  point. 

Along  the  conical  entrance  surface  neither  equation  (20)  nor 
(21)  applies,  since  the  wall  is  at  an  angle  to  the  grid  lines.  However, 
since  the  cone  surface  is  at  an  angle  of  30°  to  the  centerline,  equation 
(21)  is  used  instead  of  equation  (20).  Although  this  is  an 
approximation,  the  conical  entrance  is  the  farthest  surface  from  the 
receiving  tube  and  is  not  expected  to  affect  the  flow  in  the  inertial 
separation  region.  Therefore,  equation  (21)  is  used  along  the  conical 
wall  with  the  subscript  P  denoting  a  boundary  node  point.  The 
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subscript  NP  denotes  the  adjacent  interior  node  point  on  a  grid  line 
parallel  to  the  entrance  and  passing  through  the  point  P. 

Even  though  the  vorticity  values  along  the  conical  surface  can 
now  be  solved,  there  is  still  one  more  problem  presented  by  this 
surface.  This  problem  arises  when  the  vorticity  is  calculated  at  the 
node  points  directly  adjacent  to  the  cone.  If  the  expressions  for  AE 
and  AN  of  equation  (13)  are  examined  it  will  be  seen  that  the  stream 
function  at  point  NE  must  be  known.  We  are  assuming  in  Figure  1 
that  the  conical  surface  is  passing  through  points  E  and  N  and  point 
NE,  which  in  (1+1,  J+l),  will  lie  outside  of  the  calculation  domain. 
Therefore,  the  expression  for  AE  and  AN  must  be  revised  when  the 
vorticity  at  nodes  in  this  position  are  calculated. 

The  revision  of  AE  and  AN  can  best  be  performed  by  noting 
that  AE  and  AN  are  actually 

AE=  <('l'se-V„e)+  -  Vnc  I)  DV 


and 


(25) 


AN  =  {(vw  -  w  )+|w  -\i/  1}  DV 

lVYne  Ynw'  ,Yne  ynw  J 

The  subscripts  refer  to  the  comers  of  the  dotted  rectangle  of  Figure 

1. 

It  can  be  seen  in  Figure  1  that  the  point  denoted  by  subscript 

ne  lies  on  the  surface  of  the  cone.  Since  the  boundary  condition  for 
the  stream  function  on  the  conical  surface  is  y=0,  Vne=0  in  equations 
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(25).  If  the  other  terms  are  expressed  as  an  average  of  their  four 
surrounding  node  points,  equation  (25)  becomes: 

AE  =  {(Vg  +  Vg£  +  +  Vp)  +  lv§  +  V§e  +  Vn  +  Vp  1}  DV 

(26) 

AN  =  {(-Vjyj  -  *  V\y  ’  ¥p)  +  1‘Vjq  'Vjsjw  ‘Vy/  '  ¥p  1}  D'V 

These  equations  must  be  used  in  place  of  AE  and  AN  in  equations 
(13)  when  the  vorticity  is  calculated  at  the  node  points  one  step  from 
the  conical  surface. 

At  the  entrance,  the  flow  is  assumed  uniform  and  parallel  to 
the  impacior  axis.  Therefore,  co/r  =  0.  At  the  exit  it  is  assumed: 

/COx  /COx 

VIN=  VlNM  (27) 

On  the  centerline  both  to  and  r  are  equal  to  zero.  Hence  co/r 
becomes  indeterminate.  However,  near  the  centerline  co/r  is  an  even 
function  or  r  so  that: 

co/r  =  A  +  Br^  (28) 

At  the  centerline  r  =  0  so  (co/r)j  =  A,  where  the  subscript  1  denotes  a 


centerline  node. 
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The  values  of  A  can  be  determined  by  evaluating  equation  (28) 
at  the  two  grid  lines  adjacent  to  the  centerline.  At  the  nearest  grid 
line 


(co/r)2  s  A  +  Br22 


(29) 


and  at  the  next  grid  line 


(co/r)3  =  A  +  Brj2 


(30) 


Since  the  grid  lines  near  the  centerline  are  equally  spaced,  r2  = 
Ar  and  r3  =  2Ar  where  Ar  is  the  grid  spacing,  A  can  thus  be  solved  as: 


✓ok  4(o/r)2  -  (co/r)3 
A  *  'r'l  =  3 


(31) 


Which  is  the  iteration  equation  used  on  the  centerline. 

After  the  stream  function  has  been  determined  at  each  of  the 
node  points,  the  velocity  component  in  both  the  r  and  z  directions 
can  be  calculated. 

The  x-component  of  velocity  can  be  determined  from  the 
definition  of  stream  function: 


Vr 


=  ldJL 

~  r  dz 


(32) 


Since  the  grid  lines  are  not  equally  spaced,  node  P  is  not  equally 
distant  from  nodes  S  and  N  as  can  be  seen  in  Figure  1.  Therefore, 
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equation  (32)  is  weighted  more  in  favor  of  the  closer  grid  point  as 
follows: 


v  -  i (dA  (-p-  ZS)  I 

vr  ~  r  v9z'n  VzM  -  r  ^ 


dzJ  s 


/zN  -  zP\ 
'zN  -  zS' 


(33) 


If  the  partial  derivatives  are  put  in  finite  difference  form,  equation 
(33)  becomes: 

v  __  i_  ~  Vp>  (ZP  -  ZS>  +  (¥P  '  ¥S)  (ZN  ‘  ZP) 
r”  rp  (ZN  -  zp)  (zN  -  zS)  +  (zp  -  zS)  (zN  -  zs)  (  } 


The  z-component  of  velocity  is  derived  in  exactly  the  same  manner 
from  the  equation: 


V 


=  .1  ll 

z  r  2  r 


(35) 


to  give: 


J_  '  yp)  (rW  '  rP) 
z  ~  rp  (rE  *  rP)  (rE  *  rW ) 


(Vp  ~  ¥\y)  (rp  -  rg) 
+  (rP  •  rW  )  (rE  *  rW ) 


(36) 
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The  value  of  Vz  at  the  centerline  cannot  be  determined  from 
equation  (36)  since  rp  =  0  for  this  case.  At  the  centerline  it  can  be 

assumed  that  the  stream  function  is  an  even  function  of  r 


y  =  ar^  + 


(37) 


where  a  and  b  are  constants.  It  can  be  seen  that  b  =  vjtp  and  a  = 
VNp  *  VP 


W2 


Again  the  subscript  P  denotes  a  node  on  the  centerline 


and  NP  denotes  the  adjacent  interior  node. 

The  definition  of  the  stream  function  can  now  be  used  to  derive 
an  expression  for  Vz. 


v 

z  t  dr 


Substitute  equation  (37)  into  (38)  and  we  obtain 


(38) 


V 


z  ~ 


-  2a 


or 


2(¥NP  -  ¥P) 
z  “  rNP2 

Which  is  used  in  calculating  Vz  on  the  centerline. 


(39) 


2.3  Numerical  Technique  for  Particle  Trajectory  Calculation 


Once  rhe  flow  field  has  been  determined  as  explained  in  the 
previous  section,  the  particle  trajectories  can  be  calculated.  The 
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starting  point  for  the  particle  trajectory  calculations  is  the  values  of 
the  velocity  components  in  the  r  and  z  directions  at  each  node  point. 
Recall  that  these  velocity  components  were  the  final  result  of  the 
flow  field  calculation. 

The  governing  equation  for  the  particle  trajectory  is  the 
Newton's  law  governing  the  particle  motion: 


F=  ma 


(40) 


where  F  is  the  aerodynamic  force  on  the  particle,  m  is  the  mass  of 
the  particle  and  a  is  the  acceleration  of  the  particle.  By  assuming 
Stokes'  law  of  resistance,  the  equations  of  motion  of  a  panicle  in  the 
r  and  z  directions  are: 


2r'  3rcp.Dp  ,  dr_ 
dt'2  C  ' V  r  "  dt' 


and 


m 


3tc|aDp  ,  dz\ 
C  vvz‘3T^ 


(41) 


These  equations  can  be  made  dimensionless  by  using: 


r  = 


W 


V 


r  ~ 


(  i) 
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t' 

1  ■  w/v0 


The  equations  now  become: 


Stk  d^r  dr 

2  dt2  “  r  "dt 


Stk  d2z  dz 

2  dt2  =  Vz*dt 


where 


Stk  = 


p„cv0d2p 

9pW 


(43) 


(44) 


The  square  root  of  the  Stokes  number,  V  Stk,  is  used  as  a 
dimensionless  particle  size. 

It  should  be  noted  that  equations  (43)  are  interrelated  since 
both  fluid  velocity  components,  Vr  and  Vz  are  dependent  on  the 

position  of  the  particle,  (r,z).  Therefore,  these  two  equations  must  be 
integrated  simultaneously. 

Several  numerical  methods  are  available  for  carrying  out  this 
integration.  The  method  we  have  chosen  is  the  Runge-Kutta  fourth 
order  method,  since  it  is  self  starting  and  gives  good  accuracy. 

In  general  the  method  for  determining  the  particle  trajectory  is 
as  follows: 


1 
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1.  A  particle  of  know  Stokes  number  is  assigned  to  start  at  some 
point  at  the  virtual  impactor  inlet  with  a  velocity  equal  to  that 
of  the  fluid  at  that  point. 

2.  Linear  interpolation  formulas  are  derived  for  the  velocity 
components  Vr  and  Vz  based  on  the  values  of  the  velocity 

components  at  the  four  node  points  situated  at  the  corners  of 
the  elemental  area  in  which  the  particle  is  located. 

3.  An  estimate  is  made  of  the  particle  crossing  time,  i.e.  the  time 
taken  by  the  particle  to  cross  this  elemental  area. 

4.  Equations  (43)  are  integrated  numerically  to  determine  a  new 

particle  position  using  the  Runge-Kutta  method  with  a  time 
increment  equal  to  1  /10th  of  the  estimated  particle  crossing 
time. 

5.  A  check  is  made  to  see  if  this  new  position  of  the  particle  is  still 
within  the  same  elemental  area.  If  the  particle  is  still  in  the 
same  area,  step  4  is  repeated  with  the  same  time  increment 
until  a  grid  line  is  crossed. 

6.  When  a  grid  line  is  crossed  the  process  is  started  again  at  step 

2  for  this  new  elemental  area.  This  is  continued  until  the 

particle  reaches  the  exit  of  either  the  major  or  minor  flows  or 
strikes  a  surface. 

By  starting  particles  of  different  sizes  (Stokes  Numbers)  at 
different  points  along  the  inlet,  it  is  possible  to  construct  efficiency 
curves  for  particles  passing  through  the  receiving  tube  with  the 
minor  flow  and  out  the  side  exit  with  the  major  flow.  In  addition,  it 
is  possible  to  determine  the  percentage  of  particles  of  any  size  that 
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hit  a  surface.  These  are  considered  to  be  particle  losses  in  the 
impactor. 

2.4  Application  of  Flow  Field  and  Particle  Trajectory  Solution 

Techniques  to  Virtual  Imnactors 

The  best  application  of  the  flow  field  and  particle  trajectory 
solution  techniques  to  aerosol  transport  problems  is  on  a  case  by  case 
basis  so  the  exact  boundary  conditions  can  be  modeled.  However, 
since  virtual  impactors  are  rather  simple  instruments  with  air 
flowing  from  nozzles  into  a  receiving  tubes,  it  is  possible  to  do  a 
parametric  study  of  virtual  impactors  and  apply  the  results  from  the 
case  that  is  the  closest  to  that  of  interest.  This  was  done  by  Marple 
and  Chien  (1980)  and  the  resulting  efficiency  curves  from  their 
study  reproduced  in  figures  2  and  3.  Figure  2  shows  the  basic 
configuration  and  defines  the  parameters  while  Figure  3  shows  the 
efficiency  curves.  For  cases  which  are  different  from  those  given  in 
these  figures,  such  as  different  geometries  or  flow  rates,  the 
technique  should  be  applied  with  the  correct  boundary  conditions. 

The  efficiency  curve  shown  in  Figure  3  are  for  the  percentage 
of  particles  passing  through  the  nozzle  and  penetrating  into  the 
receiving  tube.  These  are  the  large  particle  collection  efficiency 
curves.  In  each  set  of  efficiency  curves,  one  parameter  has  been 
varied.  Therefore,  it  was  possible  to  study  the  effects  c  varying 
parameters  such  as  the  angle  of  the  inlet,  the  distance  between  the 
receiving  tube  inlet  and  the  nozzle  exit,  the  Reynolds  Number,  the 
configuration  of  the  inlet  to  the  receiving  tube  and  other 
relationships  between  the  diameters  of  the  receiving  tube  to  that  of 
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Figure  2. 
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Figure  3.  Virtual  impactor  efficiency  curves  showing  (1)  effect  of 
jet  Reynolds  number,  (b)  ratio  of  minor  to  total  flow,  (c) 
ratio  of  receiving  tube  to  jet  diameter,  (d)  jet  to  plate 
distance,  (e)  nozzle  angle  and  (0  effect  of  collection  probe 
inlet. 
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the  nozzle.  One  striking  feature  of  these  efficiency  curves  is  that  all 
parameters  except  the  minor  to  total  flow  ratio  parameter  have  only 
a  minor  influence  on  the  position  or  the  shape  of  the  efficiency  curve. 
The  reason  for  this  is  that  the  flow  from  the  nozzle  penetrates  into 
the  receiving  tube  and  all  of  the  separation  is  accomplished  within 
the  receiving  tube.  Other  parameters,  such  as  the  distance  between 
the  receiving  tube  inlet  and  nozzle  exit,  the  angle  of  the  inlet,  the 
nozzle  Reynolds  Number,  shape  of  the  receiving  tube,  all  have  little 
influence  on  the  collection  characteristics. 

The  one  parameter  that  has  a  significant  influence  on  the 
efficiency  curve  is  the  ratio  of  minor  to  total  flow.  As  the  minor  flow 
is  decreased,  the  efficiency  curve  shifts  to  the  larger  sizes  of  STK 
indicating  a  larger  cutoff  size.  In  addition,  the  cutoff  characteristics 
of  the  efficiency  curves  become  sharper.  Also  note  that  the 
efficiency  curves  asymptotically  approach  the  value  corresponding  to 
the  percentage  of  flow  passing  through  the  receiving  tube  for  the 
smaller  particle  sizes.  Therefore,  it  is  desirable  to  have  as  small  a 
minor  flow  as  possible  to  obtain  a  sharp  cutoff  characteristic. 

In  the  above  described  calculation  results,  it  is  assumed  that 
particles  impacting  on  the  receiving  tube  walls  will  stick  to  the  walls. 
Since  this  may  not  occur  in  actual  practice  when  solid  particles  are 
involved,  there  is  some  uncertainty  in  the  use  of  these  theoretical 
results  for  actual  impactor  designs.  For  these  and  other  reasons, 
impactors  designed  on  the  basis  of  theoretical  calculations  must  be 
calibrated  and  studied  experimentally. 
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3.  Design  Considerations 

The  starting  point  in  the  design  of  virtual  impactor 
concentrators  is  the  Stokes  number,  defined  as  the  ratio  of  the 
particle  stopping  distance  to  the  radius  of  the  nozzle, 

x  U 

Stk=-y  (45) 


where  Stk  is  the  Stokes  number,  z  is  the  particle  relaxation  time  and 

U  is  the  mean  jet  velocity  in  the  nozzle.  The  particle  relaxation  time, 
t,  is  a  function  of  the  particle  diameter.  Dp,  and  particle  density,  pp, 


.-PjiflL- 

“  18  n  C 


(46) 


where  C  is  the  particle  slip  correction  given  by 
C  =  1  +  1.246  Kn  +  0.42  Kn  e  '  0.87/Kn 
and  Kn  is  the  Knudsen  number  defined  by 


(47) 


(48) 


To  design  a  virtual  impactor  concentrator,  the  functional  relationship 
between  the  efficiency  of  the  concentrator  and  the  Stokes  must  be 
known.  Theoretical  and  experimental  studies  of  the  virtual  impactor 
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described  above  have  shown  that  the  collection  efficiency,  r\,  is 
primarily  a  function  of  the  Stokes  number,  St, 

T|  =  f(Stk)  (49) 

with  the  other  parameters,  such  as  the  jet  to  receiving  tube  distance, 
jet  Reynolds  number,  the  nozzle  to  receiving-tube  diameter  ratio,  the 
major  and  minor  flow  ratios,  etc.  playing  a  secondary,  though 
sometimes  important,  role.  Therefore,  the  design  of  virtual  impactor 
concentrator  can  begin  by  first  as  ming  a  constant  Stokes  number 
at  the  desired  impactor  cut-point  and  then  calculating  the 
corresponding  nozzle  diameters  and  required  jet  velocities  to  achieve 
the  desired  cut  point  point.  The  results  can  then  be  verified  by 
constructing  an  experimental  impactor  and  performing  the  necessary 
calibration  and  verification  studies.  The  design  procedure  is 
described  in  more  detail  below. 
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r=<-Mt)1/3  (51) 

for  the  nozzle  radius.  Knowing  the  nozzle  radius,  the  required  jet 
velocity,  U,  can  then  be  calculated  by  means  of  Equation  (50).  To 
determine  the  pressure  drop,  AP,  across  the  nozzle,  we  make  use  of 
the  Bernoulli's  law  and  the  definition  of  the  nozzle  discharge 
coefficient,  Cd, 


U  =  Cd 


from  which  we  obtain 


AP  = 


p  U2 
2  Cd2 


(52) 


(53) 


The  power  required  to  maintain  an  air  flow  at  the  rate  of  q  through  a 
nozzle  at  a  pressure  drop  of  AP  is 


Pth  =  q  AP 


(54) 


which  is  the  theoretical  air  power  for  operating  the  virtual  impactor 
concentrator.  The  actual  input  power,  P,  to  the  blower  or  pump  is 
higher  by  a  factor  equal  to  the  blow  or  pump  efficiency,  115, 


P  = 


Pth 

Tib 


(55) 


30 


Figure  4  shows  the  calculated  nozzle  diameter,  jet  velocity,  pressure 
drop  and  the  theoretical  air  power  for  a  single  nozzle  virtual 
impactor  concentrator  designed  for  a  cut-point  diameter  of  2.0  pm . 
The  calculation  has  been  made  by  assuming  V  St  =  0.65  and  a  nozzle 
discharge  coefficient  of  0.8.  Other  values  can  be  used  for  the 
calculation,  but  the  results  will  not  differ  substantially  from  those 
obtained. 

It  is  interesting  to  note  that  the  nozzle  diameter,  the  jet 
velocity,  the  pressure  drop  and  the  theoretical  air  power  all  increase 
substantially  with  increasing  flow  rate.  For  a  1  liters  per  minute 
(lpm)  sampling  flow  rate  a  nozzle  of  1.11  mm  in  diameter  is 
required.  The  required  nozzle  diameter  increases  to  1.11  cm  at  a  flow 
rate  of  1000  lpm.  The  corresponding  jet  velocity  increases  from 
1,740  cm/s  to  17,400  cm/s,  and  the  pressure  drop  from  1.14  to  114 
inch  H2O.  The  theoretical  air  power  incr  ases  from  4.71  mW  »o  471 

W.  More  detailed  data  on  the  design  calculations  are  given  in 
Appendix.  A. 

3.2  Multiple -nozzle  Virtual  Impactor  Concentrator 

Because  of  the  high  pressure  drop  and  power  requirement  of 
the  single  nozzle  virtual  impactor  concentrator,  practical  impactors 
with  moderate  to  high  flow  rate  requirements  generally  must  be 
designed  with  multiple  nozzles.  For  a  multiple  nozzle  impactor 
consisting  of  n  identical  nozzles  operating  in  parallel,  the  flow  rate 
through  each  nozzle  is 


Single  Nozzle  Virtual  Impactor  Concentrator 
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q  = 


ft 


(56) 


where  Q  is  total  flow  rate  through  the  multiple  nozzle  impactor.  by 
substituting  Eq.  (56)  into  Equations  (51)  and  making  use  of  Equations 
(51)  -  (54),  we  have  the  following  equations  for  the  design  of 
multiple  nozzle  impactor 

Q  x  , 

R  =  (  1/3  (57) 

n  it  o  t 


U  =  Cd 


(58) 


AP  = 


p  U2 
2  Cd2 


(59) 


Pth  =  Q  AP 


(60) 


Figures  5-8  show  the  calculated  impactor  characteristics  for  a 
virtual  impactor  concentrator  operating  with  a  flow  rate  of  10,  100, 
1000  and  22  1pm  and  a  cut  point  diameter  of  2,0  pm.  In  this 
calculation,  the  number  of  nozzles  has  been  varied  from  1  to  100. 
The  advantage  of  the  multiple  nozzle  approach  is  thus  obvious.  For 
instance,  for  a  virtual  impactor  concentrator  with  a  sampling  flow 
rate  of  22  1pm,  a  single  nozzle  design  will  result  in  a  concentrator 
having  a  pressure  drop  of  8.92  inch  H2O.  With  100  nozzles,  the 
pressure  drop  is  reduced  to  0.414  inch  H2O.  The  decrease  in  the 
theoretical  air  power  is  quite  dramatic.  A  single  nozzle  impactor 


Multiple  Nozzle  Virtual  Impactor  Concentrator  ( Q=  10 1pm ) 
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Figure  5.  Dependence  of  nozzle  diameter,  pressure  drop,  jet  velocity,  and  air  power  on  number  of 
nozzles  for  a  virtual  impactor  concentrator  (Q  =  10  1pm,  D  =  2.0  pm) 


Multiple  Nozzle  Virtual  Impactor  Concentrator  ( Q=  100 1pm  ) 
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Dependence  of  nozzie  diameter,  pressure  drop,  jet  velocity,  and  air  power  on  number  of 
nozzles  for  a  virtual  impactor  concentrator  (Q  =  100  Ipm,  D  =  2.0  urn) 


Multiple  Nozzle  Virtual  Impactor  Concentrator  (  9=  1000 1pm ) 


Dependence  of  nozzle  diameter,  pressure  drop,  jet  velocity,  and  air  power  on  number  of 
nozzles  for  a  virtual  impactor  concentrator  (Q  =  1000  Ipm.  D  =  2.0  uml 


Multiple  Nozzle  Virtual  Impactor  Concentrator  ( 9= 


Dependence  of  jet  velocity  on  number  of  nozzles  for  a  multiple  nozzle  virtual  impactor 
concentrator  (O  -  22  Ipm,  D  =  2.0  um) 
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requires  a  power  of  814  mW;  while  a  100  nozzle  impactor  requires 
only  37.8  mW  of  theoretical  air  power.  More  detailed  data  on  the 
design  calculations  are  given  in  Appendix  B. 

4.  Design  and  Evaluation  of  a  22  Liter  Per  Minute  Virtual 

Imttaqtor  -  Concentrator 

On  the  basis  of  the  above  design  considerations,  a  virtual 
impactor  having  a  sampling  flow  rate  of  22  liters  per  minute  was 
designed,  constructed  and  tested.  The  concentrator  consists  of  22 
nozzles  each  operating  at  a  flow  rate  of  1  liter  per  minute.  The 
specifications  of  the  design  are  as  follows: 


Nozzle  diameter  = 

0.048" 

Total  flow  rate  = 

22  1pm 

Major  flow  rate  = 

21  1pm 

Minor  flow  rate  = 

1.0  1pm 

Number  of  nozzles^ 

22 

A  design  drawing  for  the  concentrator  is  shown  in  Figure  9.  Detailed 
design  of  the  various  parts  of  the  system  is  given  in  Appendix  C. 

4.1  Experimental  Evaluation  and  Results 

For  the  purpose  of  the  laboratory  evaluation  of  the  virtual 
impactor  concentrator,  an  inlet  adapter  piece  is  constructed  so  that 
the  test  aerosol  can  be  introduced  into  the  concentrator  through  this 
adapter  piece.  Figure  10  shows  the  concentrator  with  this  inlet 
adapter  piece  used  during  the  test  in  place  of  the  regular  inlet. 
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The  actual  calibration  was  made  using  liquid  oleic-acid 
particles  generated  with  a  Berglund-Liu  vibrating  orifice 
monodisperse  aerosol  generator  (Berglund  and  Liu,  1973)  to 
determine  its  particle  collection  efficiency  and  loss  characteristics. 
Fluorescent  uranine-tagged  oil  drops  of  a  specific  particle  size  were 
generated  by  the  Berglund-Liu  generator,  and  were  collected  using 
the  large  particle  and  small  particle  collection  filters  shown  in  Figure 
11.  Particle  loss  inside  the  minor  flow  chamber  was  also  determined 
by  wiping  the  internal  surfaces  with  a  wet  cotton  swap  and 
measuring  the  fluorescence.  The  experimental  results  are  shown  in 
Figure  12. 

Results  of  these  experiments  show  that  the  virtual  impactor 
concentrator  has  excellent  particle  concentrating  performance 
characteristics.  The  50%  aerodynamic  was  found  to  be  3.1  pm  and  the 
internal  particle  loss  in  the  minor  flow  chamber  was  quite  small, 
ranging  from  0.5%  for  a  particle  diameter  of  2.53  pm,  to  a  maximum 
of  12.8%  at  the  particle  diameter  of  7.31  pm. 

The  pressure  drop  through  the  concentrator  was  checked  at 
various  flow  rates  and  at  the  locations  shown  in  Figure  13  using  an 
inclined  water  manometer.  The  results  are  shown  in  Figure  14.  For 
the  operating  conditions  of  21  1pm  major  flow  and  1  1pm  minor  flow, 
the  pressure  drops  are  0.9  and  0.03  in.  H2O  respectively. 

4.2  Conclusions 

The  experimental  virtual  impactor  concentrator  was  found  to 
perform  satisfactorily  with  high  collection  efficiency,  low  particle  loss 
and  low  pressure  drop.  The  experimentally  determined  cut-point 
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Figure  11.  Virtual  impactor  concentrator  with  large  particle  and 
small  particle  collection  filre  s  during  test  to  determine 
collection  efficiency  and  Ie  loss 
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Figure  12.  Collection  efficiency  and  particle  loss  characteristics  of 
virtual  impactor  concentrator 


Figure  13.  Virtual  impactor  collector  during  pressure  drop  tests, 
showing  locations  where  pressure  was  measured 
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Figure  14.  Pressure  drop  vs.  flow  characteristics  of  virtual  impactor 
concentrator 


diameter  at  3.1  pm  of  the  concentrator  was  larger  than  intended. 
However,  the  cut-point  diameter  can  be  decreased  by  increasing  the 
jet  velocity.  To  achieve  a  cutpoint  diameter  of  2.0  pm  it  is  necessary 
to  increase  the  jet  velocity  by  a  factor  equal  to  the  ratio  of  the 
relaxation  time  of  the  3.1  pm  diameter  particle  to  that  of  the  2.0  pm 
diameter.  The  ratio  is  found  to  be  2.34.  By  increasing  the  total 

flowrate  of  the  concentrator  from  22  1pm  to  51  5  1pm,  the  desired 
cutpoint  can  be  achieved.  At  this  flow  rate,  the  pressure  drop 

through  the  impactor  would  be  4.0  in.  H2O,  based  on  the  flow  vs. 

pressure  drop  measurement.  No  other  change  in  the  design  of  the 
concentrator  would  be  necessary. 

5.  V  aporizer/Pvrolvzer 

The  purpose  of  the  vaporizer/pyrolyzer  is  to  collect  airborne 
particles  from  the  concentrated  aerosol  stream  from  the  virtual 

impactor  concentrator  onto  a  suitable  substrate  and  then  vaporize  or 
pyrolyze  the  collected  particles  for  subsequent  physical  or  chemical 
analysis.  Specific  chemical  detection  techniques  that  can  be  used 

include  mass-spectroscopic  or  gas  chromatographic  techniques. 
However,  this  project  deals  with  the  vaporizer/pyrolyzer  designs 
only.  The  specific  chemical  analysis  techniques  to  be  used  are  not 
investigated. 

5.1  Bidiminaiy  Experiments 

Preliminary  experiments  were  performed  by  means  of  a  foil 
vaporizer/pyrolyzer  made  ot  nichrome  foil.  The  physical 
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characteristics  of  the  nichrome  foil  used  in  the  experiments  were  as 
follows: 


Resistivity: 

108  \iCl/cm3 

Composition: 

Chromium 

20% 

Iron 

0.5% 

Silicon 

1.0% 

Nickel 

78.5% 

Dimensions: 

0.25"  x  1.0"  x  0.001"  (Width  x  length  x  thickness) 

The  material  was  obtained  from  Hoskins  Manufacturing  Co.  (4445 
Lawton  Ave.,  Detroit,  Michigan  48208)  and  the  foil  was  mounted  on 
supporting  electrodes.  The  completed  unit  was  found  to  have  an 
electrical  resistance  of  0.345  Cl. 

The  electrical  and  thermal  characteristics  of  the  foil 
vaporizer/pyrolyzer  was  determined  by  connecting  a  low  voltage  DC 
power  supply  to  the  evaporator  and  measuring  the  electrical  current 
and  the  temperature  of  the  foil  as  the  voltage  is  varied.  The  power 
k  supply  used  has  a  voltage  range  of  0  -  36  Volts  and  current  range  of 
0-10  Amperes.  The  foil  temperature  was  measured  by  welding  a 
thermocouple  to  the  surface  of  the  foil  and  measuring  the 
thermocouple  output  directly  with  a  temperature  indicator.  The 
experimental  data  are  summarized  in  Table  1  and  the  results  are 
plotted  in  Figures  15  and  16  for  the  voltage  vs.  current  and  power 
vs.  temperature  relationships.  It  is  interesting  to  note  that  the  foil  is 
able  to  achieve  a  temperature  of  700  °C  with  an  input  power  of  28 
Watts. 
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Table  1:  Electrical  and  Thermal  Characteristics  of  Foil 

Vaporizer/Pyrolyzer 


I,  Ampere 

V.  Volt 

W,  Watt 

T,°C 

0.0 

0.00 

0.00 

25 

1.0 

0.34 

0.34 

50 

2.0 

0.69 

1.38 

94 

3.0 

1.01 

3.03 

145 

4.0 

1.38 

5.52 

205 

5.0 

1.72 

8.60 

290 

6.0 

2.07 

12.48 

398 

7.0 

2.44 

17.08 

510 

8.0 

2.76 

22.08 

615 

9.0 

3.11 

27.98 

700 

Table  2:  Vaporization  Spectrum  for  4  pL  of  0.1%  DOP 
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The  Performance  of  the  vaporizer/pyrolyzer  was  evaluated 
using  the  experimental  system  shown  in  Figure  17.  Her-;  the 
vaporizer/pyrolyzer  is  mounted  inside  a  small  glass  bottle  and 
filtered  air  at  24  or  12  1pm  is  introduced  into  the  bottle  to  flush  out 
the  vaporized  or  pyrolyzed  material.  A  TSI  3760  condensation 
nucleus  counter  (CNC)  is  used  to  monitor  the  partLics  in  the  output 
air  stream.  Since  the  CNC  is  capable  of  detecting  particles  as  small  as 
0.01  pm,  the  system  proves  to  be  extremely  sensitive  for  detecting 
the  small  amounts  of  vaporized  or  pyrolyzed  material  from  the 
vaporizer/pyrolyzer. 

Foil  Cleaning:  For  the  purpose  of  cleaning  the  foil,  the 
vaporizer/pyrolyzer  was  immersed  in  denatured  alcohol  for  about 
one  hour  and  then  dried  in  an  oven.  The  assembly  was  then 
installed  in  the  system.  Next,  23  watts  of  power  was  applied  so  that 
all  material  that  would  be  vaporized  or  pyrolyzed  below  this  power 
level  would  be  vaporized  or  pyrolyzed.  Finally,  4  pL  of  denatured 
alcohol  was  added  to  the  foil.  The  foil  was  then  heated  with  23  w  of 
power.  No  CNC  count  was  observed,  indicating  that  the  foil  was  in  fact 
clean  and  that  the  alcohol  used  was  nearly  impurity  free. 

Vaporization  Spectrum:  To  determine  the  response  of  the 
vaporizer/pyrolyzer,  a  solution  was  prepared  by  adding  0.1%  of  DOP 
(di-octyl  phthalate)  by  volume  to  the  denatured  alcohol  solvent.  4  pL 
of  this  0.1%  DOP  (di-octyl  phthalate)  solution  was  then  added  to  the 
foil.  The  electrical  current  was  then  slowly  increased.  When  the  CNC 
began  to  indicate  a  count,  the  current  increase  was  stopped  and  the 
current  was  maintained  at  that  level  until  the  CNC  count  rose  to  a 
peak  and  then  decayed  to  zero.  The  current  was  then  slowly 
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increased  and  the  above  steps  repeated  until  all  current  levels  at 
which  counts  were  observed  were  determined. 

Figure  18  shows  a  sample  spectrum  obtained  with  4  pL  of  0.1% 
DOP.  It  is  interesting  to  note  that  two  peaks  were  observed  at  power 
levels  of  2.09  W  and  3.83  W,  corresponding  to  temperature  levels  of 
approximately  110  an."!  170  °C.  At  these  temperature  levels, 

materials  in  the  DOP  were  vaporized  and  recondenses  to  form  an 
aerosol,  which  is  then  detected  by  the  CNC.  The  detailed  data  are 
given  in  Table  2. 

Flash  Vaporization:  For  the  flash  vaporization  tests,  2,  4  or  8  pL 
of  the  sample  solution  was  added  to  the  vaporizer/pyrolyzer  and  12 
W  of  power  was  applied  to  the  foil  vaporizer/pyrolyzer  suddenly. 
Power  was  maintained  at  this  level  and  aerosol  particle  counts  were 
recorded  by  the  CNC.  The  results  obtained  for  the  0.1%  DOP  are 
shown  in  Table  3.  It  is  interesting  to  note  that  the  peak 
concentration  detected  is  approximately  proportional  to  the  solution 
volume  added  to  the  foil,  indicating  a  linear  relationship  between 
vaporized  material  and  peak  aerosol  concentration  detected. 

5.2  Remits 

DOS :  Experiments  were  performed  using  both  solutions  of 
known  impurity  content  as  well  as  solutions  with  unknown  impurity 
content.  Figure  19  shows  the  spectrum  obtained  with  4  pL  of  0.1% 
DOS  (di-octyl  sebacate)  in  denatured  alcohol.  The  spectrum  is  quite 
similar  to  that  for  DOP  with  the  second  peak  occurring  at  a  power 
level  somewhat  lower  than  that  for  DOP  (3.32  W  for  DOS  vs.  3.83  W 
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Table  3:  Flash  Vaporization  for  0.1%  DOP 


Volume,  pL 

#/cc 

Run  1 

Run  2 

2 

2.2  x  104 

2.4  x  104 

4 

6.2  x  104 

6.2  x  104 

8 

1.2  x  105 

8.9  x  104 

Table  4:  Vaporization 

Spectrum  for  4  pL  of  0.1%  DOS 

#/cc 

W,  Watt  T,°C 

Run  1 

Run  2 

2.04 

3.1  x  104 

4.5  x  104 

3.32 

2.1  x  104 

3.6  x  102 

Table  5:  Flash 

Vaporization  for  0.1%  DOS 

#/cc 

Volume,  pL 

Run  1 

Run  2 

2 

3.2  x  104 

3.0  x  104 

4 

4.6  x  104 

5.3  x  104 

8 

7.0  x  105 

6.7  x  104 

51 


for  DOP).  The  data  are  given  in  Table  4.  The  results  for  the  flash 
vaporization  tests  are  given  in  Table  5. 

De-ionzed  Water:  To  test  the  impurity  in  de-ionized  water,  two 
grades  of  de-ionized  water  was  used.  One  of  the  de-ionized  water 
tested  was  known  to  be  of  high-quality  with  electrical  resistivity  in 
the  vicinity  of  18  MO.  The  second  de-ionized  water  was  considerably 
lower  in  quality,  with  much  lower,  but  unknown  electrical  resistivity. 
The  results  are  shown  in  Figure  20.  In  both  cases,  one  peak  was 
observed  at  a  power  level  in  the  vicinity  of  13  W.  The  result  for  the 
flash  vaporization  test  is  shown  in  Table  6. 

Tap  Water:  When  tap  water  from  the  laboratory  faucet  was 
used,  interesting  and  complicated  spectra  was  observed.  Figure  21 
shows  the  spectrum  for  the  tap  water.  A  total  of  six  peaks  were 
observed  corresponding  to  the  various  unknown  impurities  in  the 
tap  water.  The  results  are  given  in  more  detail  in  Table  7.  The 
corresponding  data  for  the  flash  vaporization  test  are  given  in  Table 
8. 

Finger  Print:  When  the  foil  vaporizer/pyrolyzer  surface  is 
lighted  touched  by  the  finger,  there  is  sufficient  contamination  to 
give  interesting  and  easily  observed  peaks.  The  results  are  shown  in 
Figure  22  and  in  Tables  9  and  10. 

The  above  experimental  results  are  compared  and  summarized 
in  Table  11.  The  highest  particle  count  is  obtained  with  DOP  where  4 
pL  of  0.1%  of  DOP  produced  the  highest  count  of  6  x  104  particles/cc. 
Lightly  touching  the  foil  surface  with  the  tip  of  the  finger  produced 
counts  that  are  of  the  same  order  of  magnitude.  On  the  other  hand, 
when  high  purity  DI  water  was  used,  the  counts  obtained  was 
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Table  6: 

Flash  Vaporization  of  De-ionized 

Water 

#/cc 

Run 

High  Purity 

Low  Purity 

1 

0.054 

3.732 

2 

0.193 

0.508 

3 

0.610 

0.447 

4 

0.054 

3.808 

5 

0.060 

1.040 

Mean 

0.194 

1.905 

Table  7: 

Vaporization 

Spectrum  for  Tap 

Water 

W,  Watt 

T,°C 

_ mL 

Run  1 

Run  2 

3.04 

142 

0.265 

0.803 

4.10 

176 

1.71 

1.342 

5.36 

202 

0.162 

0.521 

8.77 

300 

0.017 

0 

14.02 

440 

1.197  x  103 

1.15x  103 

19.59 

580 

1.070  x  104 

2.26x  102 

Table  8:  Flash  Vaporization  of  4  pL  of  Tap  Water 


Run 

1 

2 

_3 _ 

Mean 


#/cc _ 

2.40  x  104 
2.39  x  104 
2.82  x  104 
2.54  x  104 


Table  9:  Vaporization  Spectrum  for  Finger  Print 


Ms 


W,  Watt 

T,  °C 

Run  1 

2.67 

2.85x102 

5.07x102 

3.02 

6.77x102 

9.29xl04 

3.90 

7.24xl04 

3.38xl04 

5.28 

4.82x10* 

1.905 

17.97 

6.47x103 

9.65x102 
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Table  10:  Flash  Vaporization  of  Finger  Print 


Run 

#/cc 

1 

4.95  x 

104 

2 

5.80  x 

104 

3 

1.78  x 

105 

4 

1.44  x 

105 

Mean 

1.07  x 

105 

Table  11:  Typical  CNC  Count  (#/cc) 

(4  pL  sample  Volume,  except 

for  Various  Materials 

finger  print) 

Material 

#/cc 

0.1%  DOP 

6  x  104 

0.1%  DOS 

5  x  104 

Tap  Water 

2.5  x  104 

DI  Water  (high  purity) 

0.194 

DI  Water  (low  purity) 

1.90 

Finger  Print 

1  x  104 
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approximately  5  orders  of  magnitude  lower,  and  in  the  vicinity  of 
0.19  particles/cc. 

5.3  YaporizerVPvrolvzer  _  for  the  Virtual  Imoactor  Concentrator 

On  the  basis  of  the  experimental  work  described  above,  a 
vaporizer/pyrolyzer  was  designed  and  constructed  for  the  virtual 
impactor  concentrator  described  earlier.  The  conceptual  design  of  the 
vaporizer/pyrolyzer  is  shown  in  Figure  23  and  a  more  detailed 
design  is  shown  in  Figure  24. 

Specifically,  the  vaporizer/pyrolyzer  consists  of  a  needle-in-a- 
tube  electrostatic  precipitator  to  precipitator  the  particles  in  the 
concentrator  aerosol  stream  from  the  virtual  impactor  concentrator. 
The  particle  collecting  surface  is  in  the  form  of  a  tight  cylindrical 
helix  wound  from  nichrome  heating  wires.  The  helical  coil  of 
nichrome  wire  provide  the  collecting  surface  onto  which  particles  can 
precipitate.  After  sufficiently  amounts  of  material  has  been  collected, 
an  electrical  current  can  then  be  passed  through  the  nichrome  wire 
to  cause  the  material  to  be  vaporized  or  pyrolyzed.  The  vaporized  or 
pyrolyzed  material  can  then  be  introduced  into  a  mass  spectrometer 
for  physical  or  chemical  analysis.  This  particular  design  of  the 
vaporizer/pyrolyzer  has  the  advantage  that  the  electrostatic 
precipitator  designed  has  high  collecting  efficiency  and  the  small 
mass  of  the  nichrome  wire  results  in  short  thermal  time  constant  and 
rapid  heating  of  the  collected  sample.  In  addition,  since  particles  are 
collected  on  the  nichrome  wire  by  electrostatic  charging  and 
precipitation,  particles  would  be  attracted  to  the  collecting  surface 
and  bounce  and  re-entrainment  of  particles  are  unlikely  to  occur. 
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This  particular  design  is  specifically  recommended  for  further 
research  and  development  effort. 

6.  Summary  of  Important  Results 

The  result  of  this  study  indicates  that  sample  acquisition  for 
chemical  or  biological  detection  can  best  be  accomplished  by  the  use 
of  a  multiple  nozzle  virtual  impactor  concentrator  to  concentrate  the 
aerosol  stream  and  an  electrostatic  precipitator  to  deposit  the 
particles  onto  the  surface  of  a  nichrome  wire  to  volatilize  or  pyrolyze 
the  collected  material  for  subsequent  physical  or  chemical 
identification.  The  use  of  a  multiple  nozzle  virtual  impactor 
concentrator  will  lead  to  significant  savings  in  pressure  drop  and 
power  consumption  for  the  sample  flow  pump.  The  use  of  an 
electrostatic  precipitator  will  provide  high  collection  efficiency  and 
avoid  the  problem  of  particle  bounce.  A  sample  acquistion  system 
incorporating  the  above  mentioned  elements  has  been  designed, 
constructed  and  tested.  The  performance  characteristics  of  the 
system  have  been  found  to  be  satisfactory. 
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Figure  15.  Voltage-current  curves  for  foil  vaporizer/pyrolyzer 
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Figure  16.  Foil  temperature  as  a  function  of  input  power 


Figure  19.  Vaporization  spectrum  for  4  nL  of  0.1%  DOS  solution 
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Figure  20.  Vaporization  spectrum  for  de-ionized  water  of  high  and 
low  purity 


Figure  21.  Vaporization  spectrum  for  ordinary  tap  water 


Watt 

Figure  22.  Vaporization  spectrum  for  finger  print 


(19.95.0.241) 


Figure  24.  Design  of  electrostatic  precipitator,  vaporizer/pyrolyzer 
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Appendix  A 

Performance  of  the  Single  Nozzle  Virtual  Impactor  Concentrator 
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DP,  urn 

2 

SliD  Correction 

1.11E+00 

Relaxation  Time 

1.34E-05 

StA0.5 

6.50E-01 

Stokes  Number 

4.23E-01 

Nozzle  Coefficient 

8.00E-01 

Nozzle  Dia. 


D,  cm 


1.11E-01 


1.39E-01 


1 .59E-01 


1.75E-01 


1.89E-01 


2.01  E-01 


2.1 1  E-01 


2.21  E-01 


2.30E-01 


2.38E-01 


3.00E-01 


3.43E-01 


Jet  Veloci 


U,  cm/s 


1.74E+03 


2.19E+03 


2.50E+03 


2.76E+03 


2.97E+03 


3.16E  +  03 


3.32E+03 


3.47E+03 


3.61E  +  03 


3.74E  +  03 


4.71E+03 


5.40E+03 


4E+ 


6.40E+03 


6.80E+03 


7.16E+03 


7.48E+03 


7.78E+03 


8.06E+03 


1.02E+04 


1 .1 6E+04 


1.28E+04 


1.38E+04 


1.46E+04 


1.54E  +  04 


1.61E+04 


1.68E+04 


1.74E+04 


Pressure  Dro 


Inch  H20 


1.14E+00 


1.80E+00 


2.36E+00 


2.86E+00 


3.32E  +  00 


3.75E+00 


4.16E+00 


4.54E+00 


4.91E+00 


5.27E+00 


8.37E+00 


1.10E+01 


E+01 


1.54E+01 


1.74E+01 


1.93E+01 


2.11  E+01 


2.28E+01 


2.45E+01 


3.88E+01 


5.09E+01 


6.17E+01 


7.16E+01 


8.08E+01 


8.96E+01 


9.79E+01 


1.06E+02 


1.14E+02 


Air  Power 


Watt 


4.71E-03 


1.50E-02 


2.94E-02 


4.75E-02 


6.89E-02 


9.34E-02 


1 .21  E-01 


1 .51  E-01 


1 .84E-01 


2.19E-01 


6.94E-01 


1.37E+00 


2.20E  +  00 


3.20E+00 


4.33E  +  00 


5.60E+00 


7.00E+00 


8.52E+00 


1.02E  +  01 


3.22E+01 


6.34E+01 


1.02E+02 


1.48E+02 


2.01  E+02 


2.60E+02 


3.25E  +  02 


3.95E  +  02 


4.71  E+02 


ron 


Mt 


mti 
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Appendix  B 

Performance  of  the  Multiple  Nozzle  Virtual  Impactor  Concentrator 
with  Sampling  Flow  Rates  of  10,  100,  1000  and  22  Liters  per  Minute 
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DP,  um _ 

Slip  Correction 


Relaxation  Time 


StA0.5 


Stokes  Number 


Nozzle  Coefficient 


Q,  Ipm 


1.11  E+00 


1 .34E-05 


6.50E-01 


4.23E-01 


Flow  Rate/ 


nozzle,  Ipm 


1.00E+01 


5.00E+00 


3.33E+00 


2.50E+00 


6  1.67E+00 


1.25E+00 


1.11  E+00 


1.00E+00 


5.00E-01 


3.33E-01 


2.50E-01 


2.00E-01 


1 .67E-01 


1 .43E-01 


1.25E-01 


1.11E-01 


1.00E-01 


Nozzle  Dia 

i. 

D,  cm 

2.38 

E- 

01 

1  .89 

E- 

01 

1  .65 

E- 

01 

1 .50 

E- 

01 

1  .39 

E- 

01 

1  .31 

E- 

01 

1  .24 

E- 

01 

1  .19 

E- 

01 

1  .14 

E- 

01 

1  .1 1 

E- 

01 

8.77 

E- 

02 

7.66 

E- 

02 

6.96 

E- 

02 

6.46 

E- 

02 

6.08 

E- 

02 

5.78 

E- 

02 

5.53 

E- 

02 

5.31 

E- 

02 

5.13 

E- 

02 

,  cm/s 


3.74E+03 


2.97E+03 


2.59E+03 


2.36E  +  03 


2.1  9E+03 


2.06E+03 


1 .96E  +  03 


1 .87E+03 


1 .80E  +  03 


1 .74E+03 


1 .20E  +  03 


1.09E  +  03 


1 .02E  +  03 


9.56E+02 


9.08E  +  02 


8.68E  +  02 


8.35E+02 


8.06E+02 


Press.  Drop 

Air  Power 

inch  H20 

Watt 

5.27E+00 

2.19E-01 

3.32E+00 

1  .38  E-0 1 

2.53E+00 

1 .05  E-0 1 

2.09E+00 

8.68E-02 

7.48E-02 

1 .60E  +  00 

6.62E-02 

1 .44E+00 

5.98E-02 

1.32E+00 

1.22E+00 

5.06E-02 

1 .14E+00 

4.71  E-02 

7.16E-01 

2.97E-02 

5.46E-01 

4.51  E-01 

1  .87E-02 

3.88E-01 

1  .61  E-02 

3.44E-01 

1  .43E-02 

3.1  0E-01 

1  .29  E-02 

2.84E-01 

1 . 1  8  E-02 

2.63E-01 

1  .09  E-02 

2.45E-01 

1  .02  E-02 

#  of  Nozzle 


Multiple  Nozzle  Virtual  Impactor  Concentrator  ( Q=  10  lpm ) 


of  Nozzle 


Multiple  Nozzle  Virtual  Impactor  Concentrator  ( 9=  10 1pm ) 


of  Nozzle 


Multiple  Nozzle  Virtual  Impactor  Concentrator  ( Q=  10 1pm ) 


of  Nozzle 
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DP.  um 


Slip  Correction 


Relaxation  Time 


StA0.5 


Stokes  Number 


Nozzle  Coefficient 


Q.  lorn 


1.11E+00 


1.34E-05 


6.50E-01 


4.23E-01 


8.00E-01 


100 


Noz2le  # 


fll 


Flow  Rate/ 


nozzlejom 


1.00E+02 


5.00E+01 


3.33E+01 


2.50E+01 


2.00E+01 


1.67E+01 


1.43E+01 


1.25E+01 


1.1 1E+01 


1.00E+01 


2.50E+ 


33 


Nozzle  Dia. 

Jet  Vete>  l 

D.  cm 

5.1 3E-01 

8.06E+03 

4.07E-01 

6.40E+03 

3.56E-01 

5.59E+03 

3.23E-01 

5.08E+03 

3.00E-01 

4.71E+03 

2.82E-01 

4.44E+03 

2.68E-01 

4.21  E-t-03 

2.57E-01 

4.03E  +  03 

2.47E-01 

3.88E  +  03 

2.38E-01 

3.74E+03 

1 .89E-01 

2.97E+03 

1.65E-01 

2.59E+03 

1.50E-01 

2.36E+03 

1.39E-01 

2.19E  +  03 

1.31E-01 

2.06E+03 

1.24E-01 

1 .96E  +  03 

1.19E-01 

1.87E+03 

1.14E-01 

1.80E+03 

1.11E-01 

1.74E+03 

yress.  Dro 


inch  H20 


2.45E+01 


1.54E+01 


1.18E+01 


9.71  E+00 


8.37E+00 


7.41  E+00 


6.69E+00 


6.12E  +  00 


5.27E+00 


3.32E+00 


2.09E+00 


1.80E+00 


1.60E+00 


1.44E+00 


1 .1 4E+OOI 


Air  Power 


1.02E+01 


6.40E+00 


4.88E+00 


4.03E+00 


3.47E+00 


3.07E+00 


2.77E+00 


2.54E+00 


2.35E+00 


2.19E+00 


1.38E+00 


1.05E+00 


8.68E-01 


7.48E-01 


6.62E-01 


5.98E-01 


5.47E-01 


5.06E-01 


4.71E-01 


Multiple  Nozzle  Virtual  Impactor  Concentrator  (  9s  100 1pm ) 


of  Nozzle 


Multiple  Nozzle  Virtual  Impactor  Concentrator  ( 9=  100  lpm ) 


of  Nozzle 


Multiple  Nozzle  Virtual  Impactor  Concentrator  ( Q=  100  lpm ) 


of  Nozzle 
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DP.  um 


Slio  Correction 


Relaxation  Time 


StA0.S 


Stokes  Number 


Nozzle  Coefficient 


Q.  tom 


1.11E+00 


1.34E-05 


6.50E-01 


4.23E-01 


8.00E-01 


1000 


Noz2le  # 


Ell 


E 


Row  Rate/ 


nozzle,  tom 


1.00E+03 


5.00E+02 


3.33E+02 


ii.50E+02 


2.00E+02 


1.67E+02 


1.43E+02 


1.25E+02 


1.11E+02 


1.00E+02 


5.00E+01 


3.33E+01 


2.50E+01 


2.00E+01 


1.67E+01 


1.43E+01 


1.25E+01 


1.11  E+01 


1.00E+01 


Nozzle  Di 


D,  cm 


1.11  E+00 


8.77E-01 


7.66E-01 


6.96E-01 


6.46E-01 


6.08E-01 


S.78E-01 


5.53E-01 


5.31  E-01 


5.13E-01 


4.07E-01 


3.56E-01 


3.23E-01 


3.00E-01 


2.82E-01 


2.68E-01 


2.57E-01 


2.47E-01 


2.38E-01 


la. 

Jet  Veloci 

fiL. 

Press.  Droi 

U.  cm/ 

s 

inch  H2 

0 

1.74E+04 


1.38E+04 


1.20E+04 


1.09E+04 


1 .02E+04 


9.56E+03 


9.08E+03 


8.68E+03 


8.35E+03 


8.06E+03 


6.40E+03 


5.59E+0 


5.08E+03 


4.71  E+03 


4.44E+03 


4.21  E+03 


4.03E+03 


3.88E+03 


3.74E+03 
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1.14E+02 


7.16E+01 


5.46E+01 


4.51  E+01 


3.88E+01 


3.44E+01 


3.10E+01 


2.84E+01 


2.63E+01 


2.45E+01 


1.54E+01 


1.18E+01 


9.71  E+00 


8.37E+00 


7.41  E+00 


6.69E+00 


6.12E+00 


5.66E+00 


5.27E+00 


Air  Power 


4.71  E+02 


2.97E+02 


2.27E+02 


1.87E+02 


1.61 E+02 


1 .43E+02 


1 .29E+02 


1 .1 8E+02 


1.09E+02 


1.02E+02 


6.40E+01 


4.03E+01 


3.47E+01 


3.07E+01 


2.77E+01 


2.54E+01 


2.35E+01 


2.19E+01 


Multiple  Nozzle  Virtual  Impactor  Concentrator  ( 9=  1000 1pm ) 


of  Nozzle 


Multiple  Nozzle  VLrtural  Impactor  Concentrator  ( Q=  1000 1pm ) 


of  Nozzle 


Multiple  Nozzle  Virtual  Impactor  Concentrator  ( Q=  1000  ) 


of  Nozzle 


I 
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DP.  urn 

2 

SIId  Correction 

1.11E+00 

Relaxation  Time 

1.34E-05 

StA0.5 

6.50E-01 

Stokes  Number 

4.23E-01 

Nozzle  Coefficient 

8.00E-01 

Q,  lorn 

22 

Nozzle  # 


Flow  Rate/ 


nozzle.  tom 


2.20E+01 


1.10E+01 


7.33E+00 


5.50E+00 


4.40E+00 


3.67E+00 


3.14E+00 


2.75E+00 


2.44E  +  00 


2.20E+00 


2.00E+00 


1.83E+00 


1.69E+00 


1.57E+00 


1.47E+00 


1.38E+00 


1.29E+00 


1.22E+00 


1 .1 6E+00 


1.10E+00 


1.05E+00 


1.00E+C0 


9.57E-01 


9.17E-01 


8.80E-01 


8.46E-01 


8.15E-01 


7.86E-01 


7.59E-01 


7.33E-01 


5.50E-01 


4.40E-01 


3.67E-01 


3.14E-01 


2.75E-01 


2.44E-01 


2.20E-01 


Nozzle  Dia. 


3.10E-01 


2.46E-01 


2.15E-01 


1.95E-01 


1 .81 E-01 


1.70E-01 


1.62E-01 


1.55E-01 


1.49E-01 


1 .44E-01 


1.39E-01 


1.35E-01 


1.32E-01 


1.29E-01 


1.26E-01 


1.23E-01 


1.20E-01 


1.16E-01 


1.14  E-01 


1 .1  2E-01 


1.11 E-01 


1.09E-01 


1.07E-01 


1.06E-01 


1.05E-01 


1.03E-01 


1.02E-01 


1.01  E-01 


9.97E-02 


9.06E-02 


8.41  E-02 


7.91E-02 


7.52E-02 


7.19E-02 


6.91  E-02 


6.67E-02 


ieehjeei 


,  cm/e 


4.87E  +  03 


3.86E+03 


3.37E  +  03 


3.07E  +  03 


2.85E  +  03 


2.68E+03 


2.54E  +  03 


2.43E  +  03 


2.34E  +  03 


2.26E  +  03 


2.19E+03 


2.13E+03 


2.07E+03 


2.02E  +  03 


1.97E  +  03 


1.93E+03 


1.89E+03 


1.86E+03 


1.82E+03 


1.79E+03 


1.76E  +  03 


1 .74E  +  03 


1.71E  +  03 


1.69E+03 


1.66E+03 


1.64E  +  03 


1.62E+03 


1.60E  +  03 


1.58E+03 


1.57E+03 


1.42E  +  03 


1.32E  +  03 


1 .24E  +  03 


1.18E  +  03 


1.13E  +  03 


1.09E  +  03 


1.05E  +  03 
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Press.  Dro 


Inch  H2 


8.92E  +  00 


5.62E+0Q 


4.29E+00 


3.54E+00 


3.05E+00 


2.70E  +  00 


2.44E+00 


2.23E+00 


2.06E  +  00 


1.92E+00 


1.80E+00 


1.70E+00 


1.61E+00 


1 .54E+00 


1.47E+00 


1.40E+00 


1.35E+00 


1.30E+00 


1.25E+00 


1.21E+00 


1.17E+00 


1.14E+00 


1.10E+00 


1.07E+00 


1 .04E+00 


1.02E+00 


9.91  E-01 


9.67E-01 


9.45E-01 


9.24E-01 


7.63E-01 


6.57E-01 


5.82E-01 


5.25E-01 


4.80E-01 


4.44E-01 


4.14E-01 


Air  Power 


8.14E-01 


5.13E-01 


3.91  E-01 


3.23E-01 


2.78E-01 


2.47E-01 


2.22E-01 


2.04E-01 


1 .88E-01 


1.75E-01 


1.65E-01 


1.55E-01 


1.47E-01 


1 .40E-01 


1.34E-01 


1.28E-01 


1.23E-01 


1 .1  9E-01 


1 .14E-01 


1.10E-01 
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Appendix  C 

Design  Drawings  for  the  Virtual  Impactor  Concentrator 
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